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IONOSPHERIC MODIFICATION: 

AN INITIAL REPORT ON ARTIFICIALLY CREATED EQUATORIAL SPREAD F 


I. Introduction 


Since the discovery of a large dropout (ionospheric hole) in the 


ionospheric total electron content during the launch of NASA's Skylab 


(Mendlllo et al., 1975 a,b) much investigation has gone into the 


controlled artificial production of such holes or plasma density 


depletions in the F region ionosphere (Bernhardt et al ., 1975 


Bernhardt, 1976; Bernhardt and da Rosa, 1977; Mendillo and Forbes 


454-455, 1977, Spring Meeting Special Session on "Artificially 


Created Holes in the Ionosphere"). Typically investigators have 


contain molecular species which react with ionospheric 0 in such a 


manner that the 0 loss rate due to these reactions is several orders 


of magnitude faster than natural ionospheric loss rates for 0 


Mendillo et al ., 1975 a,b; Bernhardt , 1976). Recently, Pongratz 


et al . [1978] conducted two high explosive (88kgra mixture of nitro- 


methane and ammonium nitrate) F region chemical releases (260 and 


280 km altitude), codenamed Lagopedo, over the Hawaiian Islands 


(midlatitude) and formed a hole (^ 100 km in diameter) via the 


predicted chemical processes. 
Note: Manuscript submitted March 22, 1978. 
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Within the context of artificial hole formation much of the 
theoretical work to date (Mendillo et al . t 1975 b; Bernhardt , 1976; 
Pongratz et al ., 1978) has centered on the complicated chemical 
reactions and molecular diffusion equations for producing these holes 
(e.g., electrodynamic forces, such as polarization electric fields, 
on the ions or electrons are not considered). iu the past few years, 
experimental ( Kelly et al ., 1976; McClure et al .,1977) and numerical 
simulation studies ( Scannapieco and Ossakow , 1976) have observed 
(under appropriate conditions) the formation of natural, rising 
plasma density depletions (holes or bubbles) in the equatorial 
Spread F ionosphere. In the present paper we wish to perform a 
marriage between the two studies and examine the behavior of a large 
plasma density depletion (hole) artificially injected into the bottom- 
side nighttime equatorial F region ionosphere. 

II. Theoretical Model and Numerical Simulation Studies 

In the present simplified model we will take the approach that 
the large hole has already been created and thus represents a new 
initial condition for our previous nonlinear numerical simulation 
equatorial Spread F studies [ Scannapieco and Ossakow , 1976; Ossakow 
et al ., 1978 a,bj. In this way we do not follow the chemical kinetics 
of the actual hole making process, rather we follow the evolution of 
the hole after it has been produced. Pongratz et al . [1978] noted 
that the entire hole in Lagopedo was produced in less than 5 minutes 
and Sjolander et al . [1978] noted, using rocket in situ ion mass 
spectrometer measurements during Lagopedo, that the dominant ion in 
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the artificially created hole was 0 + . Thus, we are able to utilize 
the two dimensional (perpendicular to the magnetic field, U) numerical 
simulation code developed for natural equatorial Spread F studies in 
the collisional Rayleigh-Taylor regime, which is represented by the 
following equations (see Scannapieco and Ossakow , 197b; Ossakow et al ., 
1978 a.b) 

f? " | x *> * 7n “ -V n ‘ n o ) O) 

V-CVinnV^i) - | (g x z)-Vn (2) 

where n is the electron density (n Q refers to the background equilib¬ 
rium density; n = n Q + n^), is the polorization (induced) potential, 
v R is the usual (used in natural Spread F studies) ionospheric 
recombination rate (see Fig. 1), v ln is the usual ion-neutral 
collision frequency, g is gravity (pointing in the -y direction), 
z - B/|B1 and V refers to the x and y directions. In our simulations 
realistic altitude profiles of and v R (see Fig. 1) and n ^ (see 
Fig. 2) have been used. 

The numerical simulation plane coincided with the equatorial 
plane spanning an altitude (y) range from 172 to 452 km in one case 
and 252 to 532 km in the other case (both with Ay - 2 km) witlv an 
east-west (x) extent of 200 km (Ax - 5 km). The system of eqns. (1) 
and (2) was initialized as follows 
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n/n " 
o 


0.03 , r-r < 13 km 
o- 

exp[a(r-r o ~15)^] - 0.97 , 15 km < r-r Q < 33 km 

1.0, r-r > 35 km (3) 

o 


3 

where a - ln(1.97)/(20) . This provided a constant 972 depletion 

over a central 30 km of the hole and the depletion decreased to 

ambient density over a 20 km extent on each side (total extent of 

bubble was 70 km). The thickness of the hole at half depletion was 

•v 63 km (see for example Pongratz et al .» 1978). The results of a 

simulation with r corresponding to x * 0, y * 300 km, i.e., release 
o o o 

altitude = 300 km, are shown in Figs. 2-5, which depict contour plots 
of n/n Q (or equivalently n ]/ n 0 ) at t ■ 0, 500, 1000, and 2000 sec. 

The t *• 0 frame shows the initial state after the hole has been 
formed by chemical processes. At t * 500 sec the hole has begun to 
steepen on its topside as it rises toward the F peak (altitude * 

354 km). This corresponds to a rise velocity ^ 28 m/sec. We note 
that the initial hole is no longer circular and plasma density 
enhancements are forming to the sides of the bubble (hole). As tne 
initial depletion rises on the bottomside it displaces regions of 
higher density which fall to lower altitudes creating enhancements. 

At t * 1000 sec the top of the bubble is at the F peak and the enhance¬ 
ments are getting larger. This should be contrasted with the natural 
Spread F developing bubble, in this ionospheric profile, which took 
'v 8000 sec to reach the peak (see Scannapieco and Ossakow , 1976; 

Ossakow et al ., 1978b). At t ■ 2000 sec the hole is well past the 
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peak and haa bifurcated on lta topside. This la similar to barium 
cloud striation bifurcation phenomena (see Zabuaky et al ., 197 3; 
Oasakow et al ., 1977) and would lead to smaller scale struatures 
if the present simulation had more resolution. The hole Is now 
accompanied by large enhancements and the enhancements in turn are 
forcing depletions to lower altitudes. This picture is tar different 
from the initial circular hole at t “ 0. This figure shows a wide¬ 
spread rapid equatorial Spread F condition. 

Figures 6-9 exhibit the results for an artificial depletion 
released at 380 km altitude (x^ “ O, y^ - 180 km) in an ionosphere 
with the altitude of the F peak equal to 434 km. In tins case the 
time evolution of the hole is much more rapid. For example,between 
t ■ 0 and 50 sec the hole is rising with a velocity v 88 m/sec. The 
initial rise velocity, here and in the previous case, is consistent 
with a bubble rise velocity given by (Oasakow and Ohaturvedl, 1978) 

V B - tg/v^ n ) f ( n j/n Q )» where f (itj/n^) is an increasing function of 
the percentage depletion, n^/n^, and dependent on bubble shape. In 
particular f corresponds here to a circular-like shape. The case 
depicted in Figs. 6-9 is also to be contrasted to the natural Spread F 
case in this geometry where a bubble took ' 1000 sec to reach the 
F peak. 

III. Concluding Remarks 

We have presented numerical simulation results, in the 
collisional Rayleigh-Taylor regime, of artificial depletions injected 
into the bottomside nighttime equatorial F region. It has been shown 
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Kin. 2 Contour plot* of constant t>|/n„ nt t - t) sec (initial condition). The 
circular contours represent the artificial depletion (•) such that the outermost con¬ 
tour in a 16% depletion and the innermoat is a 96% depletion (inaide this contour 
ia a conatant 07% depletion). The larne dnahed curve represent* a plot of the am¬ 
bient electron number denaity (values on upper horizontal axis), n„, as a function 
of altitude. The vertical (y) axia represent* altitude, the lower horizontal (x) axiN 
east west ratine, and the ambient mannetic field is out of the finure (z). 








FiK. .’l - Swiii' im Ki||. (1), except t - MM) «v, The smaller dashed contours 
are contours of constant n,/n„ for plasma (tensity enhancements (♦). The 
enhancement contour represents a lift enhancement. 
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